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[i]  From  historical  hydrocasts  in  the  Okinawa  Trough 
region  of  the  East  China  Sea  (ECS),  acoustic  echo  time 
from  700  dbar  to  the  surface  shows  a  tight  relationship  with 
temperature  {T)  except  near  the  100-200-dbar  layer.  This  is 
caused  by  2nd-  or  higher-mode  baroclinic  variations. 
Significant  out-of-phase  correlation  between  upper  and 
lower  layer  T  from  the  hydrocasts  confirms  that  2nd-mode 
variations  are  strong.  Furthermore,  the  2nd-mode  variations 
are  dominant  during  the  ebb  tide  period,  which  suggests 
they  are  caused  by  2nd-mode  internal  tides  (ITs)  generated 
at  the  continental  shelf  break.  For  comparison,  we 
investigate  historical  hydrocasts  southeast  of  the  Ryukyu 
Islands,  where  no  significant  2nd-mode  variations  are 
found.  A  model  simulating  the  M2  IT  agrees  qualitatively 
with  these  observations:  it  predicts  relatively  strong  2nd- 
mode  IT  energy  in  the  ECS,  but  weak  energy  in  the 
region  southeast  of  the  Ryukyu  Islands.  Citation:  Park,  J.  H., 
M.  Andres,  P.  J.  Martin,  M.  Wimbush,  and  D.  R.  Watts  (2006), 
Second  mode  internal  tides  in  the  East  China  Sea  deduced  from 
historical  hydrocasts  and  a  model,  Geophys.  Res.  Lett.,  33,  L05602, 
doi:  1 0. 1 029/2005GL024732. 

1.  Introduction 

[2 J  The  East  China  Sea  (ECS)  shelf  slope  has  been 
predicted  to  be  the  second  strongest  M2  internal  tide  (IT) 
generation  region  among  the  world’s  continental  shelf 
slopes  [Baines,  1982],  Strong  IT  generation  in  this  region 
has  been  confirmed  by  observations  and  numerical  models 
[e.g.,  Kuroda  and  Mitsudera,  1995;  Niwa  and  Hibiya , 
2004] .  Kuroda  and  Mitsudera  [1995]  observed  the  detailed 
structure  of  ITs  in  the  ECS  using  an  underwater  sliding 
vehicle.  In  particular,  they  found  abundant  evidence  of  2nd- 
mode  IT  structures  at  the  shelf  break  in  the  central  ECS 
(29°N,  127.5°E)  during  ebb  tide  periods.  They  showed  that 
this  observation  was  consistent  with  Hibiya1  s  [1986]  model 
since  the  barotropic  tidal  current  was  faster  (slower)  than  the 
phase  speed  of  the  2nd-mode  (1  st-mode)  ITs  on  the  shelf, 
and  hence  ebb  tides  were  optimal  for  the  generation  of  2nd- 
mode  ITs  in  the  ECS. 

[\]  During  2002-2004,  an  array  of  11  pressure-sensor- 
equipped  inverted  echo  sounders  (PIESs),  was  deployed  on 
two  Kuroshio-transects  lines  (Figure  1)  to  observe  Kuroshio 
variability  in  the  ECS.  The  PIESs  measured  hourly  bottom- 
to-surface  acoustic  echo  time  (t)  and  bottom  pressure 
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during  the  two-year  deployment  period.  To  obtain  temper¬ 
ature  ( T)  and  specific  volume  anomaly  (6)  profiles  from  t 
and  sea  surface  temperature  (SST),  observed  respectively  by 
PIES  and  satellite,  we  constructed  multi-index  gravest 
empirical  mode  (MI-GEM)  fields  [Park  el  a /.,  2005]  using 
historical  hydrocasts.  However,  comparison  between  the 
MI-GEM  fields  and  the  historical  hydrocasts  reveals  that 
the  MI-GEM  fields  fail  to  capture  44%  and  52%  of  the  T 
and  8  variance,  respectively,  in  the  p  =  100-200-dbar  layer. 
This  may  be  caused  by  2nd-  or  higher-mode  variations  of 
vertical  structure,  since  t  is  insensitive  to  all  but  the  lowest 
baroclinic  mode  [Watts  and  Rossby ,  1977]. 

[4]  In  this  study,  we  demonstrate  that  2nd-mode  ITs  in  the 
ECS  are  responsible  for  failure  of  the  MI-GEM  fields  to 
capture  adequately  T and  6  variations  in  the p  =  1 00-200-dbar 
layer  as  seen  in  historical  hydrographic  measurements.  We 
verify  our  observations  with  a  numerical  simulation  of  ITs. 

2.  Data  and  Methods 

[5]  We  used  1392  historical  hydrocasts  collected  in  the 
ECS  from  1929  to  2003,  as  shown  in  Figure  1.  The  data 
were  from  the  North  Pacific  Hydrobase  [Macdonald  et  ai, 
2001]  and  the  Nagasaki  Marine  Observatory  (NMO),  Japan 
Meteorological  Agency,  and  duplicate  data  were  removed. 
We  constructed  the  MI-GEM  fields  (not  shown  here)  for  T 
and  8  as  functions  of  t70o,  sea  surface  temperature  (SST), 
and  p ,  where  t7()()  is  t  from  the  700-dbar  level  to  the 
surface.  To  test  the  MI-GEM  representation,  we  calculated 
percent  variance  captured  by  the  MI-GEM  fields  for  Tand  8 
as  a  function  of  pressure  for  each  of  these  two  variables 
[Watts  et  ai,  2001]. 

[6]  We  used  sea  level  (SL)  data  at  Naha,  Okinawa  (red 
square  in  Figure  1)  to  investigate  the  time-varying  IT 
energy.  The  SL  data  were  obtained  from  the  Japan  Ocean¬ 
ographic  Data  Center  through  their  Web  site.  Larsen  et  al. 
[1985]  showed  that  the  M2  tide  is  dominant  in  the  ECS  and 
flows  normal  to  the  continental  shelf.  From  the  M2  tidal 
chart  by  Niwa  and  Hibiya  [2004],  the  tidal  phase  difference 
between  Naha  and  the  Okinawa  Trough  region  is  within  15° 
(~0.52  h).  Since  precise  observation  times  are  given  for 
NMO  hydrocasts  from  1987  to  2003,  we  could  separate 
them  into  two  groups:  flood  tide  period  (from  2.5  hours 
after  low  tide  to  2.5  hours  after  high  tide)  and  ebb  tide 
period  (from  2.5  hours  after  high  tide  to  2.5  hours  after  low 
tide),  based  on  the  M2  tide  SL  at  Naha.  For  each  group,  we 
calculated  the  correlation  coefficients  between  120-dbar  T 
(7"i 20)  and  T  at  other  depths  to  investigate  the  possible 
relation  between  T  variation  and  the  ECS  tide.  For  compar¬ 
ison,  we  did  the  same  calculation  with  NMO  hydrocasts 
from  the  region  southeast  of  the  Ryukyu  Islands. 
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Figure  1.  The  East  China  Sea  and  adjacent  seas.  Orange 
triangles  indicate  PIES  sites  and  red  square  Naha  coastal 
tide  gauge.  Blue  and  green  dots  are  hydrocasts  with  t700  < 
0.924  s  and  >0.924  s,  respectively.  Bathymetry  contours  are 
in  meters. 

[7]  The  ocean  model  used  in  this  study  is  the  Navy 
Coastal  Ocean  Model  (NCOM)  as  described  by  Martin 
[2000].  This  model  is  similar  in  its  physics  and  numerics 
to  the  Princeton  Ocean  Model  [Blumberg  and  Mellor , 
1987],  but  uses  an  implicit  treatment  of  the  free  surface. 
The  hydrostatic  assumption  of  the  NCOM  is  not  valid  in 
small  horizontal  scales,  but  the  ITs  we  focus  on  have 
relatively  large  horizontal  scales  and  should  not  be 
significantly  affected  by  this  assumption.  For  the  numer¬ 
ical  simulations  conducted  here,  the  model  domain  cov¬ 
ered  121  —  131  °E,  23-32°N  with  horizontal  grid  resolution 
of  about  2.9  km.  The  vertical  grid  was  a-coordinate  with 
40  layers,  a  uniform  stretching  in  the  vertical,  and  a 
maximum  depth  of  5500  m.  The  initial  T  and  salinity  (5) 
were  from  Levitus  climatology,  horizontally  averaged 
over  the  model  domain.  Radiation  conditions  were  used 
at  the  open  boundaries.  The  M2  barotropic  tide  was 
forced  with  elevations  and  velocities  at  the  open  boundaries 
from  the  Oregon  State  University  global  tide  database  [Egbert 
and  Erofeeva ,  2003]  and  tidal  potential  forcing  over  the 
interior.  The  model  was  run  for  20  days  with  a  200  s  timestep. 
The  model  output  used  here  consisted  of  hourly  values 
subsampled  once  every  four  grid  points  along  the  zonal  and 
meridional  directions  during  the  last  7  days. 

[8]  We  investigated  the  spatial  distribution  of  the  1st-  and 
2nd-mode  M2  IT  energy  using  the  following  procedure  at 
each  subsampled  grid  point  of  the  IT  model.  Firstly, 
fluctuation  of  T  was  converted  into  vertical  displacement 
(Q.  Then,  complex  amplitudes  of  the  M2  harmonics  were 
calculated  for  (,  and  horizontal  baroclinic  velocity  compo¬ 
nents  (w,  v).  Secondly,  a  modal  decomposition  of  the  first  5 
vertical  modes  was  calculated  using  the  initial  T and  5,  and 
the  model  bottom  depth  (//).  Then,  the  calculated  modal 
functions  were  least-squares  fitted  to  the  M2  harmonics  of  w, 
v,  and  Q.  Finally,  using  the  first  2  modes,  we  calculated  the 
depth-integrated  1st-  and  2nd-mode  potential  and  kinetic  IT 
energy  densities  given  by  [Merrifield  and  Holloway ,  2002] 

Ep„  =  t  I  p0  [.'V2|<J2] dz,  and 

Ek„  /  p0[K|2  +  |v„|2j<fe,  (I) 


respectively,  where  n  is  the  mode  number  and  p0  the  initial 
vertical  mean  density. 

3.  Results 

[9]  Figure  2a  shows  the  percent  variance  captured  by  the 
MI-GEM  fields  for  T  and  6  as  a  function  of  pressure. 
Averaged  through  the  depth  range  from  100  to  200  dbar, 
MI-GEM  captures  only  56%  and  48%  of  the  variances  for  T 
and  6,  respectively.  Figures  2b  and  2c  show  scatter  plots  of 
t7oo  vs.  r120  and  t70o  vs.  7400,  respectively,  where  T40o  is 
temperature  at  the  400-dbar  level.  Figures  2b  and  2c  reveal 
that  t700  has  a  tight  relation  with  7400,  but  not  with  7j20, 
especially  when  t700  is  less  than  0.924  s  (blue  dots  in 
Figure  1).  This  results  in  the  low  value  of  percent  variance 
captured  around  100-200  dbar  in  Figure  2a.  The  unrepre¬ 
sented  T  and  6  variations  in  the  100-200-dbar  layer  are 
presumably  from  2nd-  or  higher-mode  baroclinic  variations, 
since  t700  is  insensitive  to  all  but  the  lowest  baroclinic- 
mode  vertical  structure,  and  the  effects  of  SST  variation  do 
not  reach  deeper  than  about  p  =  100  dbar  in  the  ECS. 

[10]  Figure  3  shows  the  correlation  coefficients  (r,  solid 
line)  between  7j20  and  Tp  with  p  between  0  and  700  dbar  for 
all  historical  hydrocasts  with  t7oo  <  0.924  s  (gray-shaded 
zone  in  Figure  2b)  in  the  ECS.  The  r  shows  a  significant 
out-of-phase  correlation  between  layers  above  and  below 
~300  dbar,  which  suggests  that  2nd  vertical-mode  variation 
is  dominant  in  this  region.  For  comparison,  Figure  3  also 
exhibits  the  vertical  structure  of  T  fluctuation,  7’  (blue 
lines),  caused  by  the  normalized  C,  of  2nd-mode  ITs  with 
a  maximum  amplitude  of  1  m  calculated  from  modal 
decomposition  using  the  mean  profile  of  historical  hydro¬ 
casts  with  t7oo  <  0.924  s  and  bottom  depths  from  700  to 
1300  dbar  at  100-dbar  intervals.  The  vertical  structures  of 
calculated  T;  are  all  consistent  with  the  observed  out-of- 


(a)  Signal  variance  captured 


Figure  2.  (a)  Signal  variance  percentage  of  total  variance 
for  T and  6  computed  from  the  MI-GEM  fields  in  the  ECS. 
(b)  Scatter  plot  of  t70o  vs.  Zj20.  Gray-shaded  zone 
indicates  highly-scattered  hydrocasts  (<0.924  s,  blue  dots 
in  Figure  1)  used  for  correlation  coefficient  calculation  in 
Figure  3.  (c)  Scatter  plot  of  t700  vs.  7400.  Thick  lines  in 
Figures  2b  and  2c  are  cubic  smoothing  spline-fitted  curves. 
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r  caused  by  normalized  2nd  mode  ^  (*C) 


Figure  3.  Correlation  coefficients  between  r12o  and  T  at 
other  pressure  levels  using  historical  hydrocasts  with  t70o  < 
0.924  s  (gray-shaded  zone  in  Figure  2b)  in  the  ECS.  Gray- 
shaded  zone  indicates  90%  confidence  intervals  for  the 
calculated  correlation  coefficients.  Blue  lines  indicate  T 
fluctuation  (T7)  caused  by  normalized  vertical  displacement 
(Q  of  2nd-mode  ITs  calculated  from  modal  decomposition 
using  the  mean  profile  of  historical  hydrocasts  (t700  < 
0.924  s)  and  bottom  depths  from  700  to  1300  dbar  at  1 00- 
dbar  intervals. 
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Figure  5.  (left)  Potential  and  (right)  kinetic  energy  density 
(JrrT2)  for  the  (top)  1  st-mode  and  (bottom)  2nd-mode  M2 
ITs  computed  using  the  model  output. 


phase  correlation  with  largest  upper-layer  T7  in  the  1 00— 
200-dbar  layer. 

[n]  Figure  4  demonstrates  how  the  2nd-mode  variations 
are  related  to  the  tides.  The  blue  and  red  lines  in  Figure  4a 
indicate  the  flood  and  ebb  tide  periods  during  a  cycle  of  the 
M2  tide  at  Naha.  Figure  4c  shows  the  r  between  T\2 o  and  Tp 
with  [)  values  from  0  to  700  dbar,  for  the  flood  (blue  line) 
and  ebb  (red  line)  tide  period  hydrocast  groups  in  the  ECS 
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Figure  4.  (a)  M2  tidal  elevation  at  Naha  coastal  tide  gauge 
station.  Red  and  blue  lines  indicate  ebb  and  flood  tide 
periods,  respectively  (see  text  for  details),  (b)  Location  map 
of  the  NMO  CTDs  in  the  ECS  (green  dots)  and  south  of  the 
Ryukyu  Islands  (cyan  dots),  (c)  Correlation  coefficients 
between  T\2o  and  T  at  other  depths  using  the  hydrocasts 
with  t 7oo  <  0.924  s  collected  in  the  ECS  by  NMO  (green 
dots  in  Figure  4b).  (d)  Same  as  Figure  4c  except  in  the 
region  south  of  the  Ryukyu  Islands  (cyan  dots  in  Figure  4b). 
Red  and  blue  lines  in  Figures  4c  and  4d  indicate  ebb  and 
flood  tide  period  hydrocasts,  respectively.  Red-  and  blue- 
shaded  zones  indicate  90%  confidence  intervals  for  the 
calculated  correlation  coefficients. 


(green  dots  in  Figure  4b).  We  use  hydrocasts  with  t70o  < 
0.924  s,  which  consist  of  323  and  368  casts  for  the  flood 
and  ebb  tide  periods,  respectively.  The  r  values  show 
stronger  2nd-mode  correlations  between  the  layers  above 
and  below  ~300  dbar  during  the  ebb  tide  period  than  the 
flood  tide  period.  For  comparison.  Figure  4d  shows  r 
calculated  in  the  region  south  of  the  Ryukyu  Islands  (cyan 
dots  in  Figure  4b).  We  again  use  hydrocasts  with  t700  < 
0.924  s,  which  consist  of  381  and  382  casts  for  the  flood 
and  ebb  tide  periods,  respectively.  No  significant  2nd-mode 
variations  are  found  in  this  region  in  either  the  flood  or  ebb 
tide  periods.  This  result  suggests  that  the  baroclinic  2nd- 
mode  variations  observed  in  the  hydrocasts  are  caused  by 
2nd-mode  ITs  generated  mainly  at  the  continental  shelf 
break.  Higher-mode  ITs  and  2nd-or-higher-mode  shorter 
time  scale  internal  waves  may  contribute  some  part  of  the 
observed  variations  in  the  100-200-dbar  layer.  Neverthe¬ 
less,  Figures  3  and  4  suggest  the  predominant  signal  not 
captured  by  t70o  arises  from  2nd-mode  ITs. 

[12]  Figure  5  shows  the  1st-  and  2nd-mode  IT  energy 
densities  calculated  from  the  NCOM  model  output.  The 
model  predicts  dominant  1  st-mode  IT  energy  as  expected, 
showing  some  “hot  spots”  of  energy  near  the  continental 
shelf  and  the  Ryukyu  Islands.  The  locations  of  hot  spots  in 
our  1  st-mode  IT  kinetic  energy  distribution  (Figure  5b) 
agree  well  with  those  in  Niwa  and  Hibiya' s  [2004]  kinetic 
IT  energy  distribution  predicted  from  a  numerical  simula¬ 
tion.  In  addition  to  the  dominant  1  st-mode  energy,  our 
model  also  predicts  strong  2nd-mode  IT  energy  in  the 
ECS.  Moreover,  although  the  model  predicts  strong  I  st- 
mode  IT  energy  in  the  region  south  of  the  Ryukyu  Islands,  it 
predicts  weak  2nd-mode  energy  there.  These  simulation 
results  agree  well  with  our  observations. 

[13]  For  detailed  comparison  between  our  observations 
and  model  output,  we  calculated  r  between  F]50  and  Tp  with 
p  values  from  0  to  700  dbar  using  T  profiles  from  the  model 
output  in  the  ECS  and  south  of  the  Ryukyu  Islands, 
indicated  by  green  and  cyan  dots  (118  and  161  grid  points) 
in  Figure  6b,  respectively.  Figure  6a  exhibits  flood  and  ebb 
tide  period  grouping  of  model  output  based  on  the  tidal 
elevation  at  Naha  in  the  model  during  an  M2  tide  cycle  in 
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Figure  6.  (a)  M2  tidal  elevation  at  Naha  from  IT  model. 
Red  and  blue  asterisks  indicate  model  output  time  for  ebb 
tide  and  flood  tide  periods,  respectively,  (b)  Location  map 
of  the  model  grid  points  in  the  ECS  (green  dots)  and  south 
of  the  Ryukyu  Islands  (cyan  dots),  where  T  profiles  are  used 
for  correlation  coefficient  calculations,  (c)  Correlation 
coefficients  between  7j 5<)  and  T  at  other  depths  in  the 
ECS  (green  dots  in  Figure  6b).  (d)  Same  as  Figure  6c  except 
in  the  region  south  of  the  Ryukyu  Islands  (cyan  dots  in 
Figure  6b).  Red  and  blue  lines  in  Figures  6c  and  6d  indicate 
ebb  and  flood  tide  period  model  output,  respectively.  Red- 
and  blue-shaded  zones  indicate  90%  confidence  intervals 
for  the  calculated  correlation  coefficients 


the  last  model  day.  For  the  flood  tide  period  group,  we  use 
hourly  T  profiles  of  model  output  from  21:00  on  model  day 
19  to  02:00  on  model  day  20.  For  the  ebb  tide  period  group, 
we  use  those  from  03:00  to  08:00  on  model  day  20.  In 
Figure  6c,  r  shows  a  significant  out-of-phase  relation 
between  layers  above  and  below  ^500  dbar  during  the 
ebb  tide  period,  but  no  such  relation  during  the  flood  tide 
period.  In  Figure  6d,  r  shows  no  significant  2nd-mode 
variations  in  either  the  flood  or  ebb  tide  periods.  These 
results  are  qualitatively  consistent  with  our  observations, 
even  though  the  depth  of  zero-crossing  and  the  r  magnitudes 
differ  from  those  in  our  observations  since  the  model  was 
am  with  a  simple  initial  stratification  independent  of  (x,  y). 


4.  Conclusion 

[14]  In  the  ECS,  the  MI-GEM  fields  fail  to  capture  44% 
and  52%  of  T and  6  variations,  respectively,  in  the  p  =  1 00 — 
200-dbar  layer.  These  unrepresented  variations  are  shown  to 
be  caused  by  strong  2nd-mode  ITs  whose  signals  appear  in 
the  hydrographic  profiles  but  not  in  the  acoustic  echo  time 
T-/00-  Using  coastal  tide  gauge  data  at  Naha,  we  have 
demonstrated  that  the  2nd-mode  ITs  are  stronger  during 
the  ebb  tide  period  than  the  flood  tide  period.  The  barotropic 
tidal  current  speed  is  close  to  the  phase  speed  of  2nd-mode 
ITs  at  the  continental  shelf  break,  optimal  for  the  generation 


of  2nd-mode  ITs  [Kuroda  and  Mitsudera,  1995].  Offshore 
tidal  flow  during  the  ebb  tide  period  can  release  the 
generated  2nd-mode  ITs  into  the  Okinawa  Trough.  We  have 
also  demonstrated  that  significant  2nd-mode  variations  are 
not  present  in  the  region  south  of  the  Ryukyu  Islands.  The 
numerical  simulation  results  of  a  3-dimensional  IT  model 
are  qualitatively  consistent  with  these  observations. 

[15]  The  principal  object  of  our  PIES  experiment  in  the 
ECS  is  to  study  the  dynamics  of  Kuroshio  variability  with 
time  periods  longer  than  a  day.  Therefore,  it  is  apparent  that 
even  though  the  MI-GEM  fields  cannot  capture  the  varia¬ 
tions  caused  by  2nd-mode  ITs,  they  will  provide  accurate 
representations  of  these  subtidal  dynamics  throughout  the 
water  column  since  the  2nd-mode  ITs  with  semi-diurnal 
period  are  responsible  for  the  large  errors  at  p  =  1 00— 
200  dbar  in  the  MI-GEM  fields. 
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